Probing unification scenarios with atomic clocks 
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We discuss the usage of measurements of the stability of nature's fundamental constants coming 
from comparisons between atomic clocks as a means to constrain coupled variations of these con- 
stants in a broad class of unification scenarios. After introducing the phenomenology of these models 
we provide updated constraints, based on a global analysis of the latest experimental results. We 
obtain null results for the proton-to-electron mass ratio = (0.68 ±5.79) x 10"" yr"^ and for the 
gyromagnetic factor ijp/gp = (— 0.72±0.89) x 10^^^ yr~^ (both of these being at the 95% confidence 
level). These results are compatible with theoretical expectations on unification scenarios, but much 
freedom exists due to the presence of a degeneracy direction in the relevant parameter space. 
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I. INTRODUCTION 

The observational evidence for the acceleration of the 
universe demonstrates that canonical theories of cosmol- 
ogy and particle physics are incomplete, if not incorrect. 
Several few-sigma hints of new physics exist, but so far 
these are smoke without a smoking gun; it's time to ac- 
tively search for the gun. 

The LHC evidence for the Higgs particle strongly 
suggests that fundamental scalar fields are among na- 
ture's building blocks. An obvious follow-up question is 
whether such fields also play a role in cosmology. They 
have been invoked to explain paradigms such as infla- 
tion, cosmological phase transitions or dynamical dark 
energy, but the most direct way to infer their presence is 
(arguably) to search for spacetime variations of nature's 
fundamental constants [l|-|3| . It is known that fundamen- 
tal couplings run with energy, and many particle physics 
and cosmology models suggest that they should also roll 
with time. One example are cosmological models with 
dynamical scalar fields, including string theory. 

Astrophysical measurements have led to claims for 0- 
0] and against 0-0 variations of the fine-structure con- 
stant a and the proton-to-electron mass ratio ^ at red- 
shifts z 1 — 3. An ongoing Large Programme at 
European Southern Observatory's Very Large Telescope 
should soon clarify matters, but a resolution may have 
to wait for a forthcoming generation of high-resolution 
ultra-stable spectrographs which include these measure- 
ments among their key science drivers. Answering this 



uestion can also shed light on the enigma of dark energy 
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Any Grand- Unified Theory predicts a specific relation 
between variations of a and /i, and therefore simultane- 
ous measurements of both provide key consistency tests. 
Our work revisits the ideas of [13, [3 and applies them in 
the same spirit as [H, , by using the most recent tests 
of the stability of fundamental constants using atomic 
clock comparisons to obtain direct constraints on the 
phenomenological parameters characterizing these unifi- 
cation scenarios. 



II. PHENOMENOLOGY OF UNIFICATION 

We wish to describe phenomenologically a class of 
models with simultaneous variations of several funda- 
mental couplings, such as the fine-structure constant 
a — /he, the proton-to-electron mass ratio ^ = mp/me 
and the proton gyromagnetic ratio Qp. The simplest way 
to do this is to relate the various changes to those of a 
particular dimensionless coupling, typically a. Then if 
a = ao(l + 5a) and 



AA: _ Aa 



(1) 



we have X = Xq{1 + kxSa) and so forth. 

The relations between the couplings will be model- 
dependent. In this section we follow [l3 . [isj . consid- 
ering a class of grand unification models in which the 
weak scale is determined by dimensional transmutation 
and further assuming that relative variation of all the 
Yukawa couplings is the same. Finally we assume that 
the variation of the couplings is driven by a dilaton-type 
scalar field (as in [l^ ) . With these assumptions one finds 
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that the variations of /i and a are related through 

^ = [0.8i?-0.3(l + 5)]— , 
fi a 



(2) 



where R and S can be taken as free phenomenological 
(mo del- dependent) parameters. Their absolute value can 
be anything from order unity to several hundreds, al- 
though physically one usually expects them to be posi- 
tive. (Nevertheless, for our present purposes they can be 
taken as free parameters to be constrained by data.) 

For our purposes it's natural to assume that particle 
masses and the QCD scale vary, while the Planck mass 
is fixed. We then have 



me 2 a 



(3) 



(since the mass is simply the product of the Higgs VEV 
and the corresponding Yukawa coupling) and 



Am„ , , ^-,AQ; 

^ = [0.8i? + 0.2(1 + 5*)] . 

mp a. 



(4) 



The latter equation is the more model-dependent one, as 
it requires modeling of the proton. At a phenomenolog- 
ical level, the choice 5* = — 1, i? = can also describe 
the limiting case where a varies but the masses don't. 
Further useful relations can be obtained jT9l - l2]| for the 
g-factors for the proton and neutron 

^ = [O.lOi?- 0.04(1 -f^)]— (5) 



^ = [0.12i?- 0.05(1 + 5)]—. (6) 

These allow us to transform any measurement of a 
combination of constants (for example a, \i and g^j into a 
constraint on the (i?, S*, a) parameter space. For atomic 
clocks, the relevant g-factors are those for Rubidium and 
Caesium, so these need to be related to those of the nu- 
cleons. The way to do this stems from [isl [l9l - l2l| . Using 
a simple shell model one has 
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-1.266 
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[0.07i?- 0.03(1 + 5)]— (7) 

a 



-0.13^+0.05(1 + 5")]—; (8) 



for our purposes in the following section, the following 
derived relation is also useful 



^jgcs/gRb) 

iacs/gRb) 



1.58 



gcs 
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A more accurate phenomenological description (moti- 
vated from experimental results and including a depen- 
dence on Qn and the spin-spin interaction) leads to 

~ [O.OWE- 0.007(1 + 5)]— (10) 
gRb a 



Clock 


VAB 


vab/vab (yr ■"") 


Ref. 


Hg-Al 


^-3.208 


(5.3 + 7.9) X 10"" 


[22] 


Cs-SFe 




(-1.9 + 0.12,ta ±2.7,j;,) X 10"" 


[23] 


Cs-H 




(3.2 + 6.3) X 10"^^ 


m 


Cs-Sr 




(1.0 + 1.8) X 10"^^ 


[25j 


Cs-Hg 




(-3.7 + 3.9) X 10~^^ 


[26] 


Cs-Yb 


1.93 

gcsfJ.a 


(0.78 + 1.40) X 10"^^ 


mi 


Cs-Rb 




(0.5 + 5.3) X 10"" 


[28] 


Cs-Yb 


1 93 


(0.49 + 0.41) X 10"" 


[29] 


Cs-Rb 




(1.39 + 0.91) X 10"^*^ 


[30] 



TABLE I: Atomic clock constraints of varying fundamental 
couplings. The second column shows the combination of cou- 
plings to which the clock is sensitive, and the third column 
shows the corresponding experimental bound. The measure- 
ments in the first seven lines were the ones used in [l^; in 
our analysis the limits from Rubidium and Ytterbium clocks 
(lines 6 and 7) have been updated to those in lines 8 and 9. 
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-0.007^+0.004(1 + 5")] 



Aa 



(11) 



Notice that these coefficients are very small, particularly 
in the last parametrization. 



III. CONSTRAINTS FROM ATOMIC CLOCKS 

By measuring the rate of two different atomic clocks 
one obtains a constraint on the relative shift of the cor- 
responding characteristic frequencies. These are propor- 
tional to certain products of fundamental couplings, and 
thus the measurement can be translated into a constraint 
of the drift of that combination. Different clock com- 
parisons are sensitive to different products of these cou- 
plings, and therefore a combined analysis can in principle 
lead to constraints on each of them. 

Table U shows the existing constraints for several pairs 
of clocks. Since the Hg-Al comparison yields a direct 
constraint on a [l^, one can use the combined dataset 
to obtain constraints in the fi-gcs plane. In [l5l [ this was 
done for the first 7 entries on the table, and we reproduce 
this (for comparison purposes) in the top panel of Fig. 
[1] The bounds coming from Rubidium and Ytterbium 
clocks (lines 6 and 7) have since improved to those in 
lines 8 and 9, and a reanalysis leads to the improved 
constraints in the bottom panel of Fig. [TJ 

Notice that the new measurements must replace the 
old ones in the analysis, as they are not independent — in 
both cases, the improved results are primarily due to a 
longer comparison time for the same set of clocks. With 
the more recent measurements the degeneracy direction 
is significantly changed. This is due to the fact that the 
Rubidium measurement (which is now the most sensitive 
one, apart from the a-only one) is not sensitive to /z. 

From this combined analysis we can also calculate the 
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FIG. 1: Atomic clock constraints on the fi-gca parameter 
space. The top panel shows the constraints obtained with 
the data discussed in [isl ]. while the bottom panel shows the 
constraints derived from the most recent data, ie using [29l.[30l| 
instead of [l^, [l^]. In both cases the one-, two- and three- 
sigma likelihood contours are plotted. Notice the change in 
the degeneracy direction. 
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FIG. 2: Atomic clock constraints on the R — S parameter 
space, using the same data as for the bottom panel of Fig. 
[1] The one-, two- and three-sigma likelihood contours are 
plotted. In the top panel the relations between the variations 
of the gyromagentic rations and a are assumed to be those of 
Eqs. (|7I8|I . while the bottom panel assumes Eqs. (|10lll|) . 



95% confidence intervals for both parameters, finding 

- = (6.8 ± 57.6) X 10~" yr"^ , (12) 



^ = (9.1 ±11.3) X 10-"yr-i; (13) 
9Cs 

the latter can be equivalently expressed in terms of gp 

^ = (-7.2±8.9) X lO-^Vr"^- (14) 
9p 

These should be compared to the result of [221] for the 
fine-structure constant (also at the 95% confidence level) 

- = (-1.7±4.9) X 10-^Vr"^ ■ (15) 
a 

There is no evidence of variations. The bound for gp is 
almost as strong the one for a, whereas the one for fj, is 
significantly weaker. This highlights the importance of 



improved experimental bounds using pairs of clocks with 
different sensitivities to a, /i and gp. 

The formalism in the previous section can be used to 
obtain constraints on the R-S parameter space, shown in 
Fig. 121 As we pointed out, the relations between the 
gyromagnetic rations and a are given by Eqs. for a 

simple shell model, while a better phenomenological de- 
scription yields Eqs. (|10mil) . Fig. [5] presents the results 
for both assumptions, quantifying the importance of this 
theoretical uncertainty: with current experimental data 
this is not critical, but as measurements improve better 
theoretical calculations will become necessary. 

Here there is a degeneracy between the two parame- 
ters, so that only a combination of them can be reason- 
ably well constrained. The degeneracy direction can be 
characterized by (S* + 1) — 2.7i? = —5 ± 15, and the al- 
lowed region has a relatively large uncertainty due to the 
fact that gp is less sensitive than to the values of R 
and S. The naively expected values (for typical GUT 
models) of i? - 30, 5 - 160 [H [3 are allowed by the 
experimental results. By separately fixing each of them 



4 



we find the following conservative bounds for the other 
i? = 61 ± 71 , assuming S = 160 (16) 

S" = 76 ± 197 , assuming i? = 30 . (17) 

These values are in agreement, at the 95% confidence, 
with both methods of calculation depicted in Fig. [2] 

IV. CONCLUSIONS 

We have considered the latest tests of the stability of 
nature's fundamental constants using atomic clocks and 
discussed their usage as a tool to constrain unification 
scenarios. A global analysis of existing measurements, 
assuming the tight bound of Rosenband et al. (22j , allows 
us to obtain separate updated constraints on and gp. 

It's worth comparing our results with the ones re- 
cently found by ^3^. Although they use a difi'erent 
parametrization that does not explicitly include Qp (and 
also a slightly different dataset), the results of both 
works agree in the case of /i: at two sigma [s^l find 
fi/n= (15±60)xl0~"'^''yr~^, while we find the marginally 
tighter = (6.8 ± 57.6) x 10"" yr'^ On the other 
hand they find a relatively weak bound for the ratio of 
the quark mass to the QCD mass scale, while we find 
a comparatively stronger bound for gp] the difference is 
due to the fact that atomic clock experiments are more 
sensitive to gp than to the ratio niq/ Kqcd- 



We carried out a first analysis of the impact of atomic 
clock measurements on the phenomenological parameters 
describing the class of varying fundamental coupling sce- 
narios under consideration: R, related to QCD physics, 
and S*, related to electroweak/Higgs physics. These mea- 
surements are only sensitive to a particular combination 
of these parameters. The experimental results are in 
agreement with theoretical expectations on unification 
scenarios. 

This R-S degeneracy can be broken by measurements 
in astrophysical systems that have different sensitivities 
to these parameters. Two examples of such systems are 
main sequence stars and neutron stars, for which parts of 
the R-S parameter space have been previously explored 
in [l^ [13] . We will discuss these issues in a future pub- 
lication. 
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